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© Stochastic Representation of Chaos Using Terminal Attractors 

Fictitious control forces stabilize what would otherwise be unstable computed trajectories. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A nonlinear version of the Liouville 
equation based on terminal attractors is 
part of a mathematical formalism for de- 
scribing postinstability motions of dy- 
namical systems characterized by expo- 
nential divergences of trajectories 
leading to chaos (including turbulence 
as a form of chaos) . The formalism can 
be applied to both conservative systems 
(e.g., multibody systems in celestial me- 
chanics) and dissipative systems (e.g., 
viscous fluids) . 

This formalism at an earlier stage of 
development was reported in “Extension 
of Liouville Formalism to Postinstability 
Dynamics” (NPO-30393), NASA Tech 
Briefs, Vol. 27, No. 9 (September 2003), 
page 56. To recapitulate: The problem is 
to predict the postinstability motions of 
a dynamic system governed by a system 
of nonlinear equations and subject to 
initial conditions. The formalism of non- 
linear dynamics does not afford means 
to discriminate between stable and un- 
stable motions: an additional stability 
analysis is necessary for such discrimina- 
tion. However, an additional stability 
analysis does not suggest any modifica- 
tions of a mathematical model that 
would enable the model to describe 
postinstability motions efficiently. The 
most important type of instability that 
necessitates a postinstability description 
is associated with positive Lyapunov ex- 
ponents. Such an instability leads to ex- 
ponential growth of small errors in ini- 
tial conditions or, equivalently, 
exponential divergence of neighboring 
trajectories. 

The development of the present for- 
malism was undertaken in an effort to 
remove positive Lyapunov exponents. 
The means chosen to accomplish this is 
coupling of the governing dynamical 
equations with the corresponding Liou- 
ville equation that describes the evolu- 
tion of the flow of error probability. The 
underlying idea is to suppress the diver- 
gences of different trajectories that cor- 
respond to different initial conditions, 
without affecting a target trajectory, 
which is one that starts with prescribed 
initial conditions. 



The Probability Distribution of Error about a target trajectory becomes flattened in uncontrolled evo- 
lution as close neighboring trajectories diverge. However, when evolution is controlled by fictitious 
stabilizing forces that create a terminal attractor in probability space, the distribution of error be- 
comes more sharply peaked about the target trajectory. 


This formalism applies to a system of 
n first-order ordinary differential equa- 
tions in n unknown dynamical (state) 
variables: 

where i is an integer between 1 and n, x, 
is one of the unknown dynamical vari- 
ables, the overdot signifies differentia- 
tion with respect to time, x is the vector 
of all the dynamical variables 
and t is time. The prescribed initial con- 
ditions are given by 



The corresponding Liouville equation 
for the evolution of the probability dis- 
tribution, P(xi,x 2 ,...x„,f), of errors in the 
initial conditions is 

|^+V-(Pf)=0 

dt 

where f is the vector of all the forcing 
functions (/i,_^,.../„). It is assumed that 
this probability distribution peaks at 
zero error (representing the prescribed 
initial conditions) . 

Fictitious control (stabilizing) forces 
[F = (P 1 ,P 2 ,...P„)] are added to the system 
of differential equations. The form of 
these forces differs from that of the ficti- 
tious stabilizing force described in the 
cited previous article: Whereas previ- 
ously, the fictitious stabilizing force was 
proportional to the gradient of the prob- 
ability density in the space of the dynam- 


ical variables, the present fictitious con- 
trol forces are functions of the differ- 
ences between expected and actual val- 
ues of the dynamical variables x,: 

1 

F i =T i[{ x i)~ x i ) 3 


where y, is a positive constant and x, is 
the expected value of x,, as given by 

(*) = f X j PdX l ...dX n 


The control forces have two impor- 
tant properties: 

• Because they vanish as x s (x), they 
do not affect the target trajectory; and 

• Because the magnitudes of their deriv- 
atives approach oo as x — > (x), they 
make the target trajectory infinitely 
stable. In other words, the target tra- 
jectory becomes a terminal attractor. 

The resulting modified system of dy- 
namical equations is j 

}c i = fi+'yi[{ x i)- x i) i 


The corresponding modified Liouville 
equation is 


f Ap 
dt ^ ■ f) '- I 


' OX: 


f,+Ti[( x i 


)~ x i) 


wherein the terminal attractors act as 
nonlinear sinks of probability. 

At the limit as x, — > (x,), one can neglect 
the real force f as being much smaller 
than the control force I'}, making it possi- 
ble to decompose the Liouville equation 
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into n independent equations and to ex- 
press Pas a product of n probabilities P t \ 



n 

i\x l ,x r ..x n )=Y\p;{x i ) 

i=X 


By use of these equations, it can be 
shown that the control forces create a 
powerful terminal attractor in probabil- 
ity space that corresponds to occur- 


rence of the target trajectory with prob- 
ability one (see figure). In configura- 
tion space (space in the sense in which 
“space” is understood in casual conver- 
sation) , the effect of the control forces 
is to suppress exponential divergence 
of close neighboring trajectories with- 
out affecting the target trajectory. As a 
result, the post-instability motion is rep- 
resented by a set of functions that de- 
scribe the evolution of such statistical 
invariants such as expectations, vari- 


ances, and higher moments of the sta- 
tistics of the state variables x f as func- 
tions of time. 

This work was done by Michail Zak of Cal- 
tech for NASA’s Jet Propulsion Labora- 
tory. Further information is contained in a 
TSP (see page 1). 

The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (818) 393-2827. 
Refer to NPO-41519. 


© Two High-Temperature Foil Journal Bearings 

These are prototypes of foil bearings for aircraft gas turbine engines. 

John H. Glenn Research Center, Cleveland, Ohio 


An enlarged, high-temperature-com- 
pliant foil bearing has been built and 
tested to demonstrate the feasibility of 
such bearings for use in aircraft gas tur- 
bine engines. At 150 mm in diameter, 
this is the largest foil bearing known to 
date. This bearing is a scaled-up version 
of a patented 100-mm-diameter foil 
bearing, augmented by coating the foil 
with a proprietary high-temperature ma- 
terial. In a companion development, a 
foil bearing as described above has been 
combined with a 150-mm-diameter ac- 
tive magnetic bearing to make a hybrid 
foil magnetic bearing. 

Foil bearings are attractive for use in 
some machines in which (1) speeds of 
rotation, temperatures, or both exceed 
maximum allowable values for rolling- 
element bearings; (2) conventional lu- 
bricants decompose at high operating 
temperatures; and/or (3) it is neces- 
sary or desirable not to rely on conven- 
tional lubrication systems. In a foil 
bearing, the lubricant is the working 
fluid (e.g., air or a mixture of combus- 
tion gases) in the space between the 
journal and the shaft in the machine 
in which the bearing is installed. At no 
or low speed, the shaft is supported at 
by a spring-loaded foil journal lining. 
Once the shaft is rotating rapidly 
enough, the hydrodynamic and viscous 
forces exerted by the flow of working 
fluid between the foil and the shaft 
force the foil away from the shaft, so 
that the shaft becomes supported by a 
film of the working fluid. 

The present enlarged, high-tempera- 
ture foil bearing has been tested at 
speeds up to 27,000 rpm (at 150 mm di- 
ameter, corresponding to a surface 
speed of 212 m/s) and at temperatures 


in excess of 1,200 °F (>649 °C). These 
speed and temperature limits exceed 
those of rolling-element bearings by sev- 
eral fold. 

The hybrid foil magnetic bearing was 
conceived to take advantage of the 
strengths of the foil and the active mag- 
netic bearing while utilizing each bear- 
ing to compensate for the weakness of 
the other, for the overall purpose of ob- 
taining high load capacity at all speeds 
and temperatures (see figure). The ac- 
tive magnetic bearing exhibits excellent 
performance at low speed, where the 
surface coating on the foil bearing has 
limited load capacity. The foil bearing 


exhibits excellent performance at high 
speed, where the active magnetic bear- 
ing can fail in response to shocks and 
other transient disturbances. 

Unlike a conventional active mag- 
netic bearing, the hybrid foil magnetic 
bearing can operate without need for a 
separate protective auxiliary/backup 
bearing. In case of failure of the active 
magnetic bearing in the hybrid foil 
magnetic bearing, the foil bearing 
plays the role of the backup bearing, so 
that a rotor can continue to run on the 
foil bearing alone and then come down 
to a safe stop. The hybrid foil magnetic 
bearing exhibits both the high load ca- 
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